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We present numerical solutions to the extended Doering- 
Constantin variational principle for upper bounds on the en- 
ergy dissipation rate in plane Couette flow, bridging the en- 
tire range from low to asymptotically high Reynolds num- 
bers. Our variational bound exhibits structure, namely a 
pronounced minimum at intermediate Reynolds numbers, and 
recovers the Busse bound in the asymptotic regime. The most 
notable feature is a bifurcation of the minimizing wavenum- 
bers, giving rise to simple scaling of the optimized variational 
parameters, and of the upper bound, with the Reynolds num- 
ber. 

PACS numbers: 47.27.Nz, 03.40. Gc, 47.20.Ft 



Rigorous results are rare in the theory of turbulence. 
There are no analytical solutions to the Navier-Stokes 
equations for fully developed turbulence, and numerical 
simulations can only cope with flows at comparatively 
low Reynolds numbers, due to the enormous number of 
degrees of freedom involved. The recent formulation by 
Doering and Constantin of a novel variational principle 
for computing rigorous upper bounds on quantities char- 
acterizing turbulent flows has therefore met with 
considerable interest. When resorting to a variational 
principle, one does no longer aim at solving the equa- 
tions of motion exactly, but rather uses these equations 
to derive inequalities that bound the relevant quantities, 
such as the rate of energy dissipation or other transport 
properties. The hope is to obtain variational inequalities 
that are, on the one hand, still technically manageable, 
but capture essential features of the full dynamical fluid 
system on the other. 

This spirit has been very much alive already some 25 
years ago in the Howard-Busse theory |^|-^|. The re- 
newed excitement about such methods stems from an 
obvious question: Does the new principle provide new 
physical insight? 

In this Letter, we answer this question for energy dis- 
sipation in plane Couette flow. We have devised a nu- 
merical scheme that allows us to exhaust the Doering- 
Constantin variational principle, with the extension in- 
troduced in Ref. \7^, for the entire range from low to 
asymptotically high Reynolds numbers. Intriguingly, we 
find that the solution to this intricate problem is or- 
ganized by a simple feature: the minimizing wavenum- 
bers bifurcate at a certain intermediate Reynolds num- 
ber, thereby determining the asymptotic scaling of the 
optimal upper bound on the energy dissipation rate. As 
a consequence of this bifurcation, this upper bound ex- 



hibits some structure at about those Reynolds numbers 
where typical laboratory shear flows become turbulent. 
We start from the equations of motion 



dtu + u • Vu + Vp = v Au, 
V-u = 



(1) 
(2) 



for an incompressible fluid confined between two in- 
finitely extended, rigid plates; v is the kinematic viscosity 
and p the kinematic pressure. The lower plate, coinciding 
with the plane z = of a Cartesian coordinate system, is 
at rest, whereas the upper one at z = h is sheared with 
constant velocity U in positive x-direction. This yields 
the no-slip boundary conditions 



u(x, y, 0, t) = 0, u(x,y,h,t) = Ux 



(3) 



(x is the unit vector in a;-direction) ; periodic boundary 
conditions are imposed in x- and y-direction. The time- 
averaged rate of energy per mass dissipated in the peri- 
odicity volume f2 is given by 




( d 3 u rf 



(4) 



Our aim is to compute a rigorous upper bound on the 
long-time limit of the non-dimensionalized dissipation 
rate 



c £ (Re) 



lim 



(5) 



where Re = Uhly is the Reynolds number. 

Let us recall that below the energy stability limit 
ReES ~ 82.65 we have the exact identity c £ (Re ) = Re -1 . 
Moreover, the laminar flow gives rise to the lower bound 
c E (Re) = Re -1 for all Re fm. In order to calculate an 
upper bound on c e for Re > ReES, we employ the back- 
ground flow approach pioneered by Doering and Con- 
stantin. Decomposing the actual velocity field u into a 
stationary, divergence-free "background flow" U (not to 
be confused with the physical mean flow), which has to 
carry the physical boundary conditions, and the devia- 
tions v from that flow, 



u(x,f) = U(x) + v(x,t) 



(6) 



inserting this decomposition into Eq. (Q), and utilizing 
the equations of motion (|l|)-(||) when eliminating the de- 
viations v, one derives an inequality that bounds c e in 



1 



terms of U | 
of the form 



We confine ourselves to background flows 



u(x) = [/#:)*, 



(7) 



where £ = z/h, and introduce a dimensionless balance 
parameter a > i, which weights the contribution to et 
that stems from the cross-term containing both back- 
ground flow and deviations. This parameter constitutes 
an additional degree of freedom and thus improves the 
bound. The bounding inequality then translates into a 
variational principle for a and the profiles (j)(C) @ : 



c £ (Re ) < min 

V ' ~ <A,o>l 



1 



4(o-i; 



D{4>] 



Re 



where D{4>} denotes the profile functional 



i 

£{0}ee |dC[0'(C)] 2 -l. 



(8) 



(9) 



The essential point now is that the candidate profiles <p 
are restricted by a spectral constraint: for each admissi- 
ble <j>, all eigenvalues A of the linear eigenvalue problem 



AV 



-2/i 2 AV 




VP, 



= v-v, 



(10) 



have to be positive; the eigenvectors V(x) have to sat- 
isfy the same homogeneous boundary conditions as the 
deviations v(x, t). The number R introduced here is the 
rescaled Reynolds number, 



R 



1 



Re. 



(11) 



Computationally, the main task is to calculate, for ev- 
ery single candidate profile </>, that rescaled Reynolds 
number R c {<fi} for which the lowest eigenvalue of jlO| ) 
passes through zero. When a is adjusted such that the 
r.h.s. of the inequality (g) is minimized for fixed <fi and 
Re, subject to the condition ([ll]), this <f> produces via 
(|J) an upper bound on c e for all Re < R c {4>}. That is, 
R c {(f>} marks the highest Re up to which <f> is admissible 
to the variational principle Exploiting the periodic 
boundary conditions, each eigenvalue A is labeled by two 
wavenumbers, k x and k y . We first keep these wavenum- 
bers fixed and compute that number R = Ro{4>}(k x , k y ) 
for which the lowest eigenvalue Afc^/c becomes zero, and 
then determine R c {4>} by minimizing over all k x and k y . 
The details of this tedious procedure will be given in 
Rcf. 0. 

We use variational profiles of the type displayed in 
Fig. [l]: boundary layer segments of width 5, modeled by 
polynomials of order n, are connected by a linear piece 
with slope p, such that the profiles are n—1 times contin- 
uously differentiable at the matching points z/h = 5 and 



z/h = 1 — 5. Thus, we have the three variational param- 
eters (5, n, and p. These profiles have not been chosen 
ad libitum, but have been shown to correctly capture the 
asymptotics of the dissipation bound in a closely related 
model problem |n[] . 

The optimized profiles resulting from the variational 
principle are depicted in Fig. Q Considering the response 
of the optimized profile parameters to the increase of Re, 
we can clearly distinguish five different regimes, (i) Up 
to Re = ReES we have the laminar regime with </>(£) = C- 
(ii) For ReES < Pe < Re\ ss 160 the laminar profile is 
deformed. The parameters 6 and n remain at their values 
0.5 and 3, respectively, while the slope p decreases from 
1 to almost 0. (iii) In the following regime, Re\ < Re < 
Re2 ~ 670, boundary layers develop. Here n still remains 
fixed, p increases again, while S decreases to its minimal 
value 0.14. (iv) Then, for Re 2 < Re < Re 3 w 1845, 
n increases dramatically from 3 to 34, thus steepening 
the profile in the immediate vicinity of each boundary, 
thereby effectively generating a new internal boundary 
layer within each boundary segment. As a consequence, 
S increases back to 0.5, so that the boundary segments 
finally join together, (v) For Re > Re% the variational 
parameters obey simple scaling laws: S = 0.5 remains 
fixed, the profile slope at the boundary becomes 



0'(O) - n oc Re, 
while the slope at the midpoint is given by 



/,'(!) ^pocRe^ 1 . 



(12) 



(13) 



The size of the internal boundary layers, i.e., the range 
close to the boundaries where the optimized profiles are 
steep, is measured by l/<f>'(0) and hence proportional to 
Re" 1 . It needs to be emphasized, however, that the op- 
timal variational profiles do not necessarily resemble the 
physically realized mean flow profiles. 

A key for understanding the above scaling is provided 
by the behavior of the minimizing wavenumbers. Re- 
markably, the minimizing k x is zero for all Re, reflecting 
the fact that there is no distinguished length scale in 
shear direction. Even more important, the minimizing 
k y bifurcates at 



Re, 



460, 



(14) 



as shown in Fig. 0. This means that for Re > Res the 



two lowest eigenvalues of Eqs. (10) are degenerate and 
pass through zero simultaneous y; it is a key property 
of the variational principle that both corresponding min- 
ima of i?o{0}(O, k y ) adopt identical values. The lower k v - 
branch {k Vt i) approaches a constant value for Re — * oo, 
namely, the very same k y that also determines the energy 
stability limit, whereas the upper branch (k y ^) scales 
proportionally to Re. Both findings are intimately con- 
nected to the power-law behavior of the profile expressed 
by Eqs. (|l2|) and (|l3|). For Reynolds numbers above Res 
two different length scales, k~\ and k~\, appear on the 
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stage, which can be identified [[10| as the extension of the 
profile's flat part in the interior and the effective widths 
of the (internal) boundary layers. The regime (iii) can 
thus be considered as a cross-over regime from deformed 
profiles with no definite length scales to profiles with well- 
developed boundary layers. 

The optimal variational bound on the energy dissipa- 
tion rate c e that results from these ingredients is depicted 
in Fig. ||. Immediately above ReES this upper bound 
cj" increases, exhibits a maximum, and then decreases 
roughly proportional to Re^ 1 ^ 4 in regime (iii). The oc- 
currence of the bifurcation changes this behavior: the 
bound passes through a distinct minimum at those Re 
where the internal boundary layers start to develop, and 
finally ascends to its asymptotic value 

lim c-(JZe) = 0.01087(1). (15) 

Re — >oo 

This asymptotic bound has to be compared to the bound 
calculated by Busse in the framework of his Optimum 
Theory [|5],[| , which refers to the limit Re — ► oo and reads 

lim c E {Re) < 0.010. (16) 

Re — >oo 

The prediction that the bound furnished by the varia- 
tional principle (|8|) should coincide in the limit Re — ► oo 
with the Busse bound ([!(]) has been made by Ker- 
swell |l2|| , who could bring the variational principle into 
a form to which Busse's so-called multi-a solutions can 
be applied. Thus, Kerswell's conclusion had to rely on all 
assumptions inherent in the Optimum Theory, whereas 
we have actually constructed a rigorous solution to the 
variational principle (||). In this way, we could not only 
confirm the correctness of Busse's asymptotic result, but 
provide a rigorous bound on c £ for all Re. Our bound 
shows non-trivial structure: a pronounced minimum fol- 
lowed by a i?e-range between 1000 and 1800 in which the 
bound's curvature changes its sign. This occurs in regime 
(iv), i.e., at about those Re where one observes the onset 
of turbulence in typical laboratory shear flows. 

In order to illustrate the importance of the proper im- 
plementation of the spectral constraint, we compare in 
Fig. H the bound obtained in this work to the bound de- 
rived previously by Doering and Constantin |l|]l3| with 
the help of elementary functional estimates. These es- 
timates are convenient to evaluate, but over-satisfy the 
spectral constraint and result in the structureless bound 
c E (Re) < l/(8\/2) « 0.088 for Re > 8-y/2. The im- 
provements due to both the introduction of the balance 
parameter a and the evaluation of the less restrictive ac- 
tual spectral constraint amounts to almost an order of 
magnitude, namely to a factor of 8 in the asymptotic 
regime. All upper bounds in Fig. || become asymptoti- 
cally independent of the Reynolds number, which means 
that they are in accordance with the scaling implied by 
classical turbulence theories. 

To conclude, we have found a scheme for solving the 
extended Doering-Constantin variational principle [00] 



for upper bounds on energy dissipation in plane Cou- 
ette flow. The bound established here is rigorous for 
all Reynolds numbers, recovers the Busse bound in 
the asymptotic regime, and represents in the important 
regime of intermediate Re the best upper bound that has 
been calculated so far. Moreover, we have pinned down 
the feature that organizes the solutions to the variational 
principle, namely the bifurcation of the minimizing wave- 
numbers illustrated in Fig. ||. Based on these results, it 
is now possible to obtain similar bounds of comparable 
quality also for other flows of interest. 

However, comparing the bound with the experimental 
data indicated in Fig. ||, it is clear that this is not the 
end of the story. The variational bound lies an order of 
magnitude above the dissipation rates measured by Rei- 
chardt and Lathrop et al. [ p^[ , and these data do 
not seem to approach a constant value, but still decrease 
with Re. If one could improve the variational principle 
such that this decrease is captured, and if it did persist 
asymptotically, one could make definite statements about 
possible intermittency corrections to classical scaling [[l6) . 
To proceed along these lines is possibly one of the greatest 
challenges in the rigorous theory of turbulence. 
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FIG. 1. Variational test profiles (j>(C) versus £ = z/h. 



FIG. 2. Metamorphosis of the optimized variational pro- 
files with the increase of the Reynolds number. We have de- 
picted the most important i?e-range on a logarithmic scale, 
beginning with the energy stability limit ReES and ending in 
the scaling regime Re > Re^. 



FIG. 3. Minimizing wavenumber(s) k y belonging to the up- 
per bound on c E displayed in Fig. ^[ The wavenumbers of the 
upper branch correspond to the inverse boundary layer thick- 
ness of the profiles, while those of the lower branch reflect the 
inverse extension of the flat interior profile part. 



FIG. 4. Bounds on c e for the plane Couette flow. Points 
denote the upper bound lZ{Re ) derived from the variational 
principle (^); the solid line on the left is the lower bound 
Cg(iie) = Re" 1 . The asymptotic value of the upper bound, 
liniRe->oo c^(Re ) = 0.01087(1), lies slightly above, but within 
the uncertainty span of Busse's asymptotic result ( |l6|) . 



FIG. 5. A synopsis of bounds on c e and experimental 
data. - Slanted straight line on the left: lower bound 
Ce_(Re) — i?e _1 . — Topmost horizontal line: upper bound 
obtained by Doering and Constantin in Refs. jj],[l3| with the 
help of an over-restrictive profile constraint and piecewise lin- 
ear profiles; c^(Re) » 0.088 for Re > 11.32. — Heavy dots: 
upper bound obtained in this work from the variational prin- 
ciple (^) supplemented by the actual spectral constraint (p"o[), 
cf. Fig. @; c7(7?e) -> 0.01087(1). - Joining dashed line: 
asymptotic upper bound ( fl6| ) derived by Busse in Refs. 
c7(i?e) — > 0.010(1). The shaded area denotes the estimated 
uncertainty of this bound. — Triangles: experimental dis- 
sipation rates for the plane Couette flow measured by Rei- 
chardt j^j. — Circles: experimental dissipation rates for the 
Taylor-Couette system with small gap as measured by Lath- 
rop, Fineberg and Swinney [ flo| . 
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